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Analysis  of  steady- state  and  transient  photoconductivity  measurements  at  room  temperature 
performed  on  c-axis  oriented  GaN  nanowires  yielded  estimates  of  free  carrier  concentration,  drift 
mobility,  surface  band  bending,  and  surface  capture  coefficient  for  electrons.  Samples  grown 
(unintentionally  ft-type)  by  nitrogen-plasma-assisted  molecular  beam  epitaxy  primarily  from  two 
separate  growth  runs  were  examined.  The  results  revealed  carrier  concentration  in  the  range  of 
(3-6)  X  1016  cm  3  for  one  growth  run,  roughly  5  X  1014 —  1  X  1015  cm  3  for  the  second,  and  drift 
mobility  in  the  range  of  500-700  cm2/(V  s)  for  both.  Nanowires  were  dispersed  onto  insulating 
substrates  and  contacted  forming  single- wire,  two-terminal  structures  with  typical  electrode  gaps  of 
—3-5  jmm.  When  biased  at  1  V  bias  and  illuminated  at  360  nm  (3.6  mW/cm2)  the  thinner  (^100 
nm  diameter)  nanowires  with  the  higher  background  doping  showed  an  abrupt  increase  in 
photocurrent  from  5  pA  (noise  level)  to  0.1-1  julA.  Under  the  same  conditions,  thicker  (151-320 
nm)  nanowires  showed  roughly  ten  times  more  photocurrent,  with  dark  currents  ranging  from  2  nA 
to  1  julA.  With  the  light  blocked,  the  dark  current  was  restored  in  a  few  minutes  for  the  thinner 
samples  and  an  hour  or  more  for  the  thicker  ones.  The  samples  with  lower  carrier  concentration 
showed  similar  trends.  Excitation  in  the  360-550  nm  range  produced  substantially  weaker 
photocurrent  with  comparable  decay  rates.  Nano  wire  photoconductivity  arises  from  a  reduction  in 
the  depletion  layer  via  photogenerated  holes  drifting  to  the  surface  and  compensating  ionized 
surface  acceptors.  Simulations  yielded  (dark)  surface  band  bending  in  the  vicinity  of  0.2-0. 3  V  and 
capture  coefficient  in  the  range  of  10_23-10-19  cm2.  Atomic  layer  deposition  (ALD)  was  used  to 
conformally  deposit  ^10  nm  of  A1203  on  several  devices.  Photoconductivity,  persistent 
photoconductivity,  and  subgap  photoconductivity  of  the  coated  nanowires  were  increased  in  ah 
cases.  TaN  ALD  coatings  showed  a  reduced  effect  compared  to  the  A1203  coated  samples. 

©  2010  American  Institute  of  Physics,  [doi:  10. 1063/1 .3275888] 


I.  INTRODUCTION 

Wide-bandgap  Ill-nitride  nano  wires  (NWs)  are  gaining 
increased  attention  for  applications  involving  nanophotonics 
and  nanoelectronics.  GaN  NW  structures  may  be  grown 
by  a  variety  of  methods  that  are  now  briefly  summarized. 
Laser  catalyzed  growth  has  been  used  to  grow  [1120] 
(a-axis)  oriented  NWs.1  Core-shell  GaN  NW  heterostructure 
light-emitting  diodes  (LEDs)  have  been  realized  where  cata¬ 
lyst  grown  a- axis  cores  were  subsequently  overgrown  with 
InGaN/AlGaN  shells.2  However,  extensive  transmission 
electron  microscopy  (TEM)  analysis  on  a- axis  GaN  NWs 
reveals  that  they  can  suffer  from  a  high  density  of  extended 
defects  and  basal  plane  stacking  faults.4  By  contrast,  [0001] 
(c-axis)  oriented  GaN  NWs  grown  by  nitrogen-plasma- 
assisted  molecular  beam  epitaxy  (PAMBE)  have  been  shown 
to  be  essentially  free  of  strain  and  defects  when  examined  by 
x-ray  diffraction,  TEM,  photoluminescence  (PL),  and  nano¬ 
mechanical  studies. 5-9  The  characteristically  low  strain  and 
low  defect  density  of  these  PAMBE  grown  c-axis  NWs  have 
been  observed  earlier  by  other  workers  (Ref.  10  and  refer- 
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ences  therein),  and  dense  arrays  of  such  structures  are  being 
used  to  realize  a  new  class  of  GaN  LED  devices.11,12 

An  important  feature  of  GaN  NWs,  given  a  particular 
doping  level  and  surface  band  bending,  is  that  they  can  be 
nearly  or  fully  depleted  of  charge  carriers  and  may  conduct 
in  a  space-charge  limited  fashion  in  certain  circumstances.13 
Indeed,  it  is  the  control  of  the  depletion  layer  thickness  that 
enables  held  effect  transistor  (FET)  sensors  based  on  these 
structures  to  function  via  electrostatic  gating  by  surface 
adsorbates.14,15  Control  of  the  depletion  layer  thickness  is 
also  important  for  the  operation  of  these  NW  devices  as  pho- 
toconductive  optical  detectors  as  will  become  evident  within 
the  context  of  this  paper.  Back-gated  NW  FET  devices  are 
commonly  used  to  estimate  transport  properties  of  the  NWs 
themselves  (e.g.,  carrier  concentration  and  mobility).3  How¬ 
ever,  in  our  experience  such  approaches  for  estimating  the 
transport  properties  of  NWs  can  be  unreliable  since  we  find 
these  devices  often  display  gate  hysteresis  and  drift.16  It  is 
thus  desirable  to  establish  alternative  approaches  for  estimat¬ 
ing  transport  properties  in  NW  structures;  in  this  regard,  pho¬ 
toconductivity  (PC)  can  offer  new  insights. 

©  2010  American  Institute  of  Physics 


0021  -8979/201 0/1 07(3)/03431 8/1 4/$30.00 


107,  034318-1 


Downloaded  12  Feb  2010  to  132.163.53.130.  Redistribution  subject  to  AIP  license  or  copyright;  see  http://jap.aip.org/jap/copyright.jsp 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

NOV  2009  2' REPORT  TYPE 

3.  DATES  COVERED 

00-00-2009  to  00-00-2009 

4.  TITLE  AND  SUBTITLE 

Steady-state  and  transient  photoconductivity  in  c-axis  GaN  nanowires 
grown  by  nitrogen-plasma-assisted  molecular  beam  epitaxy 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

University  of  Colorado, Department  of  Chemistry  and 

Biochemistry, Boulder, CO, 80309 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS (ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

ARSTRATT 

1 8 .  NUMBER  1 9a.  NAME  OF 

OF  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Same  aS 

unclassified  unclassified  unclassified  Report  (SAR) 

14 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


034318-2  Sanford  et  al. 


J.  Appl.  Phys.  107,  034318  (2010) 


PC  studies  of  GaN  NWs  have  been  reported  by  a  number 
of  authors.  Calarco  et  al.13  described  dark  conductivity  and 
PC  studies  as  a  function  of  NW  diameter  for  c-axis  oriented 
unintentionally  doped  (UID)  n- type  NWs  grown  by  PAMBE. 
They  concluded  that  the  surface  band  bending  (in  the  dark) 
was  approximately  0.55  V  for  samples  thick  enough  to  be 
undepleted.  They  showed  that  the  PC  mechanism  arose  from 
photoexcited  holes  drifting  to  the  NW  surface,  screening  the 
depletion  field,  thereby  increasing  the  diameter  of  the  neutral 
region  and  increasing  the  NW  conductivity.  It  is  interesting 
to  note  that  the  basic  mechanism  of  photoexcited  holes 
screening  a  depletion  field  to  promote  gain  in  planar  III- 
nitride  photoconductive  detectors  has  been  discussed 
earlier.17 

PC  studies  of  deep  levels  in  metal-organic  chemical- 
vapor  deposition  (MOCVD)  Ni-catalyst  grown  UID  a- axis 
GaN  NWs  were  performed  by  Armstrong  et  al.  where  it 
was  argued  that  majority  carrier  photoemission  from  deep 
levels  within  the  depletion  region  resulted  in  increased  NW 
conductivity.  PC  studies  of  Si-doped  ft-type  PAMBE  grown 
c-axis  NWs  were  reported  by  Richter  et  al.19  where  they 
observed  that  Si  incorporation  would  change  the  NW  growth 
morphology  from  essentially  uniform  in  thickness  to  tapered 
and  conical.  They  noted  that  fitting  dark  current  data  on  a 
variety  of  NWs  should  in  principle  permit  determination  of 
carrier  mobility.  Mansfield  et  al  also  reported  dark  conduc¬ 
tivity  and  PC  studies  in  PAMBE  grown  Si-doped  c-axis  GaN 
NWs,  compared  the  results  with  back-gated  FET  behavior  on 
the  same  samples,  and  numerical  simulations  of  FET  behav¬ 
ior.  Moreover,  Mansfield  et  al.  estimated  the  uncertainty  in¬ 
volved  when  attempting  to  reconcile  free  carrier  concentra¬ 
tion,  drift  mobility,  and  band  bending  using  dark  current  and 
PC  data,  and  also  showed  that  the  incorporation  of  Si  in  the 
PAMBE  growth  process  did  not  appear  to  substantially  alter 
the  NW  growth  morphology,  in  contrast  to  the  preceding 
reference. 


II.  EXPERIMENTAL  PROCEDURE  AND  DATA 
ANALYSIS 

A  number  of  parameters  and  conventions  appearing  in 
this  paper  are  defined  as  follows:  the  background  ft-type  do¬ 
nor  density  is  Nd  and  it  is  assumed  that  such  donors  are  fully 
ionized,  nL  is  the  integrated  free  charge  per  unit  length  along 
a  NW,  nph  is  the  density  of  photogenerated  holes,  D{  is  the 
diameter  of  a  cylinder  that  inscribes  a  NW  of  hexagonal 
cross  section  with  Ri=Di/ 2,  Dc  and  Rc  are  the  same  respec¬ 
tive  parameters  for  a  circumscribing  cylinder,  Dm  is  the  di¬ 
ameter  of  a  hypothetical  NW  of  circular  cross  section  with 
Rm=Dm/2,  R{  is  the  radius  of  a  cylinder  with  an  equivalent 
cross  sectional  area  as  a  hexagon  defined  by  the  circumscrib¬ 
ing  cylinder  with  radius  Rc ,  the  power  intensity  illuminating 
a  sample  at  a  wavelength  X  (in  nanometers)  is  P(X),  the  total 
steady- state  current  (dark  current  and  photocurrent)  under 
constant  illumination  is  Ip,  the  initial  transient  photocurrent 
at  the  instant  illumination  ceases  is  I°p,  the  decaying  current 
that  persists  after  illumination  ceases  is  Ip(t),  the  maximum 
“flatband”  current  conducted  in  a  NW  with  zero  depletion 


layer  thickness  is  If ,  and  the  dark  current  that  remains  after 
persistent  photoconductivity  (PPC)  has  relaxed  is  Id.  Other 
parameters  are  defined  where  they  appear. 

This  portion  of  the  paper  is  organized  as  follows:  Section 
II  A  describes  NW  sample  growth  and  fabrication.  Section 
II B  presents  an  experimental  study  of  steady-state  dark  cur¬ 
rent  and  photocurrent  and  focuses  on  samples  from  growth 
run  B738.  Section  II  C  presents  a  cylindrical  depletion  model 
which  is  applied  to  samples  displaying  nonzero  dark  current. 
Section  II  D  presents  a  model  describing  the  time  decay  of 
persistent  PC.  Section  II  E  illustrates  the  effect  on  PC  that 
result  when  NWs  are  coated  with  atomic  layer  deposition 
(ALD)  films.  Subbandgap  PC  is  also  described  in  this  sec¬ 
tion. 

A.  NW  growth  and  processing  to  fabricate  two- 
terminal  photoconductive  test  structures 

The  c-axis  oriented  GaN  NWs  used  in  this  study  were 
grown  by  PAMBE  on  Si(lll)  substrates.  Details  on  the 
growth,  spectroscopic  studies,  and  mechanical  properties  of 
NWs  grown  in  our  laboratory  are  given  elsewhere. 5-9  No¬ 
table  in  these  background  studies  is  that  the  PAMBE  grown 
wires  are  essentially  free  of  strain  and  defects.  Indeed,  TEM 
reveals  that  basal  plane  stacking  faults  are  only  found  in  the 
interfacial  region  where  the  NWs  nucleate  on  the  Si  substrate 
and  are  not  found  elsewhere  in  the  NWs.  Depending  on  the 
precise  growth  conditions,  Dc  for  the  wires  examined  ranged 
from  roughly  100  to  700  nm.  The  NWs  produced  in  a  par¬ 
ticular  growth  run  were  generally  well  formed  elongated 
hexagonal  cylinders  that  were  roughly  uniform  in  diameter. 
The  growth  rate  was  in  the  range  from  100  to  300  nm/h  and 
NWs  were  grown  to  lengths  ranging  from  5  to  20  jam. 

NWs  were  ultrasonically  stripped  off  of  the  Si(lll) 
growth  substrates  in  an  isopropanol  bath  and  dispersed  with 
a  micropipette  onto  fused  quartz  substrates  and  oxidized  Si 
substrates.  Photoconductive  test  structures  were  fabricated  in 
two  separate  formats:  (i)  wires  lying  flat  on  the  dispersal 
substrates  and  electrically  contacted  at  their  ends  and  (ii) 
wires  electrically  contacted  at  their  ends  but  suspended  hori¬ 
zontally  above  the  substrates  by  a  few  hundred  nanometers, 
thus  forming  NW  bridge  structures.  Examples  of  contacted 
NWs  dispersed  flat,  and  forming  a  bridge  structure,  are 
shown  in  Figs.  1(a)  and  1(b),  respectively.  The  NWs  are 
quite  flexible  and,  as  illustrated  in  Fig.  1(c),  in  certain  in¬ 
stances,  a  bridging  wire  would  bend  down  into  contact  with 
the  substrate  without  breaking  or  losing  electrical  continuity. 
Various  test  structures  fabricated  from  NWs  taken  from  the 
separate  PAMBE  growth  runs  are  summarized  in  Table  I.  For 
NWs  of  comparable  Dc  taken  from  a  common  growth  run, 
there  was  no  significant  difference  in  photoconductive  be¬ 
havior  for  test  structures  composed  of  wires  lying  flat  on 
dispersal  substrates  or  in  the  form  of  bridge  structures.  The 
fabrication  scheme  for  NW  electrical  contacts  is  now  de¬ 
scribed. 

Contacts  were  formed  by  electron  beam  deposition  of  20 
nm  of  Ti  followed  by  200  nm  of  Al.  Subsequent  lift-off  pho¬ 
tolithography  yielded  a  dense  array  of  indexed  electrodes 
with  single  NWs  occasionally  bridging  an  electrode  pair.  The 
samples  were  then  processed  by  rapid  thermal  annealing 
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(C) 


FIG.  1 .  (a)  Single  NW  lying  flat  on  a  fused  quartz  substrate  and  electrically 
contacted  at  the  two  ends.  This  image  corresponds  to  sample  7,  as  enumer¬ 
ated  in  Table  I.  (b)  Example  of  a  NW  bridge  structure.  The  sample  rests 
roughly  200  nm  above  the  surface  of  the  fused  quartz  substrate,  (c)  NW 
flexibility  is  illustrated  in  this  image.  During  dielectrophoresis  alignment, 
the  portion  of  the  NW  between  the  pads  sagged  and  adhered  to  the  substrate 
via  van  der  Waals  forces.  Despite  appearance  to  the  contrary,  this  device  was 
unbroken  and  exhibited  electrical  continuity. 

(RTA)  in  a  flowing  atmosphere  [1.4  SLM  (SLM  denotes 
standard  liters/minute)]  of  5%  H2/95%  Ar.  The  RTA  tem¬ 
perature  schedule  was  ramp  from  ambient  to  100  °C  in  5 
min,  dwell  at  100  °C  for  5  min,  ramp  to  500  °C  in  10  s, 
dwell  at  500  °  C  for  60  s,  and  cool  to  ambient.  The  assembly 
of  NW  bridge  structures  employed  methods  of  ac 
dielectrophoresis21  to  orient  the  NWs  between  Ti  electrode 
pads  that  were  photolithographically  defined  on  the  sub¬ 
strates.  In  these  experiments,  NWs  in  a  solvent  suspension 
were  dispersed  onto  the  substrate  while  an  ac  voltage  was 
applied  between  the  electrode  pads.  As  the  solvent  evapo¬ 
rated,  the  wires  would  come  to  rest  suspended  between  ad¬ 
jacent  electrode  pads,  thus  forming  bridges  with  the  wires 
elevated  a  few  hundred  nanometers  above  the  substrate.  Con¬ 


tacts  to  ends  of  these  bridging  NWs  were  subsequently  made 
using  another  lithography  step  and  the  same  Ti/Al  metalliza¬ 
tion  and  RTA  procedure  described  previously.  The  electrical 
conductivity  enabled  by  our  metallization  and  annealing 
scheme  is  expected  to  arise  from  the  diffusion  of  N  from  the 
NW  into  the  Ti  layer,  thereby  producing  a  TiN  layer  with  a 
reduced  work  function  compared  to  the  pure  metal.  Addition¬ 
ally,  the  N-depleted  region  in  the  GaN  under  the  metal  con¬ 
tact  is  expected  to  have  enhanced  ft-type  conductivity  via  the 
formation  of  N  vacancies.22  With  this  processing  scheme,  a 
conservative  estimate  of  the  resulting  specific  contact 
resistivity  is  in  the  range  of  (1  - 10)  X  10-5  II  cm2.  Anticipat¬ 
ing  the  discussion  that  follows  below,  we  estimated  that  Nd  is 
approximately  (3-6)  X  1016  cm-3  for  growth  batch  B738 
and  approximately  5  X  1014—  1  X  1015  cm-3  for  batch  B992. 
Following  Mohney  et  al.25  and  assuming  an  estimated  range 
for  drift  mobility  /z  of  500-700  cm2/(V  s),  we  find  that  the 
contact  resistance  is  typically  ten  times  less  than  the  NW 
resistance  considering  the  typical  ranges  of  NW  diameters 
and  electrode  gaps  given  in  Table  I. 

We  note  that  our  estimate  for  /z  is  consistent  with  the 
best  report  (that  we  are  aware  of)  for  MBE  grown  GaN  films. 
Nikishin  et  al 26  reported  a  Hall  mobility  of 
800  ±  100  cm2/(V  s)  for  crack-free  molecular  beam  epitaxy 
(MBE)  grown  GaN  with  a  n-type  background  of  (2-3) 
X  1016  cm3.  Using  a  conservative  estimate  for  the  Hall  fac¬ 
tor  of  roughly  1.2-1. 3  for  GaN,27  our  presumed  range  in  drift 
mobility  /z  is  reasonable  and  is  also  consistent  with  the  esti¬ 
mated  results  that  follow  below.  Taken  together,  these  obser¬ 
vations  and  estimates  lead  us  to  conclude  that  most  of  the 
voltage  is  dropped  across  the  NW,  and  the  resistance  of  the 
contacts  may  be  ignored  in  the  subsequent  analysis. 


TABLE  I.  Summary  of  samples  used  in  the  present  study.  All  relevant  parameters  are  defined  in  the  text,  with 
the  exception  that  the  “A”  and  T”  designations  refer  to  ALD-coated  samples  with  A1203  and  TaN,  respectively. 
For  samples  15  and  16,  two  NWs  share  common  electrode  gaps  as  indicated. 


Sample 

Dc 

(nm) 

h 

(H 

h 

b 

v. 

(V) 

Iald  (A,  T) 

Batch/die 

1 

320 

3.4 

1.1  jjlA 

5.4  julA 

0.5 

12.2  fiA  (A) 

B738/F3  1-7-1 

2 

311 

4.0 

1.1  jjlA 

11.4  julA 

1 

19.1  /i A  (A) 

B738/F3  1-7-5 

3 

224 

3.9 

54  nA 

1.6  fiA 

0.3 

B738/T157  1-4-5 

4 

214 

4.3 

94  nA 

4.7  fiA 

1 

8.6  /tA  (A) 

B738/F3  1-10-1 

5 

151 

4.2 

2.2  nA 

1.7  julA 

1 

B738/T155  II-3-6 

6 

105 

3.5 

Noise 

0.34  julA 

1 

1.1  fiA  (A) 

B738/F3  1-6-1 

7 

108 

3.8 

Noise 

0.29  julA 

1 

B738/P6  II-8-3 

8 

113 

4.2 

Noise 

0.34  julA 

-1 

2.2  fiA  (A) 

B738/F3  1-9-3 

9 

680 

3.1 

1.1  nA 

0.13  julA 

-4 

2.4  fiA  (A) 

B992/P8  1-8-4 

10 

147 

7.3 

Noise 

30.6  nA 

-4 

34.6  nA(A) 

B992/P8  II-9-3 

11 

160 

4.8 

Noise 

1.0  nA 

-4 

10.0  nA(T) 

B992/F8  II- 1-7 

12 

111 

5.8 

Noise 

1.0  nA 

4 

53.4  nA(A) 

B992/P8  III-4-1 

13 

531 

3.4 

Noise 

7.0  nA 

-4 

11.0  nA(T) 

B992/F8  II-2-1 

14 

583 

3.7 

Noise 

0.12  julA 

4 

2.2  /iA  (T) 

B992/F8  1-6-7 

15 

160 

3.4 

Noise 

0.33  julA 

-1 

0.87  /*A(T) 

B738/P5  1-1-1 

76 

16 

200 

4.5 

Noise 

16.6  nA 

1 

10.0  nA(T) 

B992/F8  1-2-1 

150 
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FIG.  2.  Schematic  of  PC  apparatus.  The  major  components  are  as  follows: 
LM:  lamp/monochromator  assembly,  LI:  collimating  lens,  M:  turning  mir¬ 
ror,  L2:  focusing  lens,  PI:  electrical  probe,  P2:  electrical  probe,  S:  sample 
under  test,  YA:  programmable  current/voltage  supply,  and  C:  control 
computer. 


B.  Dark  current  and  steady-state  photocurrent  in  GaN 
NWs 

The  NW  test  structures  were  examined  with  a  computer- 
automated  current- voltage  (/-V)  probe  apparatus  that  permit¬ 
ted  simultaneous  illumination  in  the  visible-UV  band  while 
the  I-V  characteristics  were  being  scanned  and  recorded.  All 
measurements  were  performed  at  room  temperature  in  ambi¬ 
ent  air.  A  schematic  of  the  apparatus  is  illustrated  in  Fig.  2. 
The  light  source  was  a  75  W  Xe  arc  lamp  dispersed  by  a  150 
mm  focal  length  monochromator  that  operated  with  a  spec¬ 
tral  bandwidth  of  roughly  5  nm.  The  monochromator  output 
was  focused  using  fused  silica  optics  to  a  maximum  illumi¬ 
nation  intensity  of  P(360)  =  3.6  mW  cm-2  for  near/above- 
bandgap  excitation  of  a  sample  under  test.  Higher  UV  inten¬ 
sities  were  provided  by  a  HeCd  laser  operating  at  325  nm 
which  provided  a  maximum  illumination  intensity  P(325) 
=  475  W  cm-2  with  the  optics  available. 

Starting  first  with  a  discussion  of  the  thicker  NWs  that 
displayed  nonzero  dark  current  (Table  I,  samples  1-5,  growth 
run  B738),  examples  of  the  dark  and  steady-state  PC  behav¬ 
ior  are  now  presented.  A  field  emission  scanning  electron 
microscopy  (FESEM)  image  of  sample  1  is  given  in  Fig.  3(a) 
and  Id-V  (dark)  and  Ip-V  [maximum  P(360)  illumination] 
data  for  this  sample  are  illustrated  in  Fig.  3(b).  Graphs  of  Ip 
as  a  function  of  P(360)  and  P( 325)  for  samples  3  and  5  are 
given  in  Figs.  3(c)  and  3(d),  respectively.  The  power  inten¬ 
sity  was  varied  by  means  of  neutral  density  filters  and  the 
bias  was  0.1  V  in  both  cases.  The  bias  was  chosen  to  keep 
the  maximum  Ip  from  exceeding  a  few  microamperes  since 
devices  would  often  fail  at  higher  currents.  The  data  illus¬ 
trated  in  Fig.  3(c)  indicate  that  Ip  rises  steeply  with  increas¬ 
ing  illumination  intensity  until  P(325)~10  W  cm-2,  but 
rises  with  significantly  reduced  slope  thereafter.  This  effect  is 
interpreted  as  a  screening  of  the  depletion  field  by  photoge¬ 
nerated  holes  as  follows. 

The  NW  geometry  offers  a  confinement  structure  that 
allows  one  to  equate  time  resolved  photoluminescence 
(TRPF)  lifetime  to  minority  carrier  lifetime  r,  as  described 
by  Ahrenkiel.  We  have  separately  performed  TRPF  on 
NWs  of  comparable  dimensions  from  these  growth  batches 
and  typically  find  r~  1  ns  in  ambient  air.  The  density  of 
photogenerated  holes  nph  under  constant  illumination  inten¬ 
sity  is  given  approximately  by  nph~  P(\)  a(\)  r/  ho).  The  ab¬ 
sorption  coefficient  at  325  nm  is  given  by  a( 325) —  1.2 
X105  cm"1.29  For  P(325)  =  100  W/cm2,  nph« 2 
X  1016  cm-3.  The  data  illustrated  in  Fig.  3(c)  suggest  that  for 
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FIG.  3.  (a)  FESEM  image  of  sample  1.  (b)  Ip-V  and  Id-V  for  sample  1  under 
UV  illumination  (3.6  mW/cm2  at  360  nm)  and  in  the  dark,  respectively,  (c) 
Combined  graphs  of  7p-P(360)  and  Ip-P(325)  for  sample  3  under  a  constant 
bias  of  0.1  V.  Illumination  is  at  360  nm  for  intensities  <3.6  mW/cm2  and  at 
325  nm  for  greater  intensity  levels.  The  same  data  are  plotted  on  both  linear 
and  semilog  scales  in  order  to  illustrate  the  transition  to  the  regime  where 
the  depletion  field  is  substantially  screened.  This  is  estimated  to  occur  for 
P(325)  -  100  W/cm2.  (d)  Combined  graphs  of  Ip-P(36 0)  and  Ip-P(325)  for 
sample  5  under  a  constant  bias  of  0. 1  V.  Illumination  scheme  and  simulta¬ 
neous  linear  and  semilog  representation  are  the  same  as  for  (c).  The  transi¬ 
tion  to  the  regime  where  the  depletion  field  is  substantially  screened  is  also 
estimated  to  occur  for  P(325)~  100  W/cm2. 

P(325)  ~  100  W/cm2,  nph  becomes  large  enough  to  substan¬ 
tially  screen  the  depletion  field  and  allow  the  neutral  con¬ 
ducting  region  to  expand  resulting  in  increased  NW  conduc¬ 
tivity.  This  approximate  analysis  also  reveals  that  the  surface 
density  of  photogenerated  holes  available  to  screen  the 
depletion  field  is  ~2X  1011  cm-2  at  this  illumination  inten¬ 
sity.  As  discussed  in  more  detail  below,  for  the  samples  con¬ 
sidered  in  this  study,  a  photogenerated  surface  hole  density  in 
the  range  of  (2-8)  X  1011  cm-2  should  fully  screen  the 
depletion  field  so  that  Ip  ~  If.  However,  a  photoexcited  hole 
density  of  roughly  2  X  10 10  cm-2  should  partially  screen  the 
depletion  field  such  that  Ip  will  fall  within  90%-95%  of  If. 
Hence,  we  interpret  the  inflection  illustrated  in  Fig.  3(c) 
where  P(325)~l-10  W/cm2  as  a  transitional  regime 
whereby  for  P( 325)  <  1  W/cm2,  PC  is  dominated  by  partial 
screening  of  the  depletion  field  by  photogenerated  holes 
while  for  P( 325)  >10  W/cm2,  the  structure  behaves  more 
as  a  conventional  photoconductor  in  the  sense  that  the  den¬ 
sity  of  photogenerated  electrons  becomes  comparable  or  ex¬ 
ceeds  Nd  with  increasing  P( 325).  Similar  behavior  is  seen  for 
sample  5,  as  illustrated  in  Fig.  3(d),  although  this  sample 
displays  comparatively  low  dark  current  and  is  nearly  de¬ 
pleted  in  the  dark. 

The  representative  behavior  of  the  thinner  NWs  (de¬ 
pleted  in  the  dark)  from  batch  B738  is  illustrated  in  Fig.  4(a) 
for  sample  7.  This  figure  illustrates  Ip-V  data  for  various 
constant  values  of  P(360).  Figure  4(b)  illustrates  Ip  for  this 
sample  under  a  constant  bias  of  0.8  V  as  a  function  of  vary- 
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FIG.  4.  (a)  Ip-V  for  sample  7  at  a  wavelength  of  360  nm  taken  at  various 
values  of  illumination  intensity.  The  magnitude  of  the  current  is  plotted  for 
negative  voltage  swings.  A:  0.65  julW  cm-2,  B:  4.3  julW  cm-2,  C: 
41.0  gW  cm-2,  D:  0.35  mW  cnT2,  and  E:  3.6  mW  cnT2.  (b)  Ip-P( 360) 
and  Ip-P(325)  for  sample  7  under  a  constant  applied  bias  of  0.8  V.  Experi¬ 
mental  limitations  prevented  attaining  the  maximum  value  for  7J(325) 
=475  W/cm2  that  was  illustrated  in  Figs.  3(c)  and  3(d). 


ing  levels  of  P(360)  and  P( 325)  up  to  a  maximum  of  roughly 
P( 325)  ~10W  cm-2.  Experimental  limitations  precluded  il¬ 
luminating  this  sample  at  the  maximum  value  of  P( 325) 
=475  W/cm2  as  described  for  samples  3  and  5.  Sample  6  is 
an  example  of  a  bridge  structure  formed  from  a  thinner  NW. 
A  FESEM  image  of  this  sample  is  shown  in  Fig.  5(a).  The 
I-V  characteristics  under  UV  illumination  and  in  the  dark  are 
shown  in  Fig.  5(b).  The  results  illustrating  the  behavior  of 
samples  6  and  7  are  generally  observed  for  the  thinner  NWs 
(whether  they  are  bridge  structures  or  lying  flat  on  sub¬ 
strates),  e.g.,  IpHd~  104-105.  Comparable  photoconductive 
behavior  in  MBE  grown  c-axis  GaN  NWs  depleted  in  the 
dark  was  also  given  by  Calarco  et  al P 

The  logarithmic  scale  used  to  display  the  magnitude  of 
Ip  in  Figs.  4(a)  and  5(b)  partially  obscures  the  fact  that  the 
Ip-V  plots  for  the  thinner  NWs  of  this  batch  (samples  6-8) 
are  sometimes  asymmetric  and  not  perfectly  linear,  whereas 
both  the  Ip-V  and  Id-V  behaviors  of  the  thicker  samples 
(samples  1-5)  are  usually  perfectly  linear.  Precedent  for  this 
trend  is  also  seen  in  the  results  of  Ref.  13  for  PAMBE  grown 
GaN  NWs  that  apparently  contain  a  n-type  background  car¬ 
rier  concentration  that  is  comparable  to  what  we  observe.  We 
argue  that  this  behavior  is  not  unexpected  in  the  thinner 
samples  for  the  following  reasons.  Given  the  contact  fabri¬ 
cation  scheme,  the  portion  of  the  NW  beneath  and  in  the 
immediate  vicinity  of  the  metal  contact  should  be  quite  con¬ 
ductive  as  previously  discussed.  On  the  other  hand,  the  por¬ 
tion  of  the  wire  extending  away  from  the  metal  contact  is 
fully  depleted  in  the  dark  and  partially  depleted  under  UV 
excitation.  Hence,  a  band  offset  should  exist  between  the 


FIG.  5.  (a)  FESEM  image  of  NW  bridge  (sample  6).  The  bridging  NW  is 
partially  obscured  by  an  adjacent  nonbridging  NW  as  indicated,  (b).  Ip-V 
and  Id-V  for  sample  6.  UV  illumination:  3.6  mW  cm-2  at  360  nm. 


FIG.  6.  Schematic  of  conduction  and  valence  band  bending  for  partially 
depleted  cylindrical  NW  in  the  dark  at  equilibrium.  The  diameter  Dm 
=  2 Rm.  The  n-type  NW  contains  a  fully  ionized  donor  density  Nd.  The  sur¬ 
face  supports  a  density  Ns  of  ionized  acceptors  that  trap  conduction  elec¬ 
trons,  resulting  in  a  cylindrically  symmetric  depletion  region  and  upward 
bending  bands.  The  Fermi  level  is  “pinned”  near  midgap  at  the  surface,  as 
shown.  The  neutral  region  of  the  NW  of  diameter  2 r0  that  will  conduct 
electrical  current  is  indicated.  The  insulating  depletion  region  supports  the 
upward  band  bending  shown.  The  maximum  band  bending  of  the  conduc¬ 
tion  band  between  the  undepleted  region  and  the  NW  surface  is  defined  as 
4>m.  The  surface  potential  is  defined  by  <t>s=Ec-Ef  evaluated  at  the  surface, 
as  shown. 


more  highly  conductive  GaN  directly  beneath  the  metal  con¬ 
tact  and  the  axially  adjacent  portion  of  the  NW  protruding 
away  from  the  contact.  Variations  in  contact  processing  con¬ 
ditions,  slight  tapering  of  the  wires,  etc.,  would  lead  to  un¬ 
equal  barriers  at  either  end  of  the  NW.  We  speculate  that  the 
apparent  partial  Schottky-like  behavior  observed  for  some  of 
the  thinner  NWs  under  UV  illumination  and  the  variation  in 
I-V  asymmetry  with  UV  illumination  intensity  as  seen  in 
Figs.  4(a)  and  5(b)  may  arise  from  these  effects. 


C.  Analysis  of  steady-state  dark  current  and 
photocurrent  in  GaN  NWs 

1.  Cylindrical  model  to  simulate  depletion  in  NWs  of 
hexagonal  cross  section 


As  illustrated  in  Fig.  6,  we  assume  that  the  NW  may  be 
approximated  by  a  cylinder  of  diameter  Dm  with  uniform 
ft-type  doping  density  Nd  leading  to  a  free  carrier  concentra¬ 
tion  of  the  same  density.  The  NW  is  partially  or  fully  de¬ 
pleted  by  the  presence  of  surface  acceptor  states  of  density 
Ns,  resulting  in  an  undepleted  conducting  neutral  region  of 
radius  rQ  within  the  wire.  The  maximum  upward  bending  of 
the  conduction  band  between  the  undepleted  neutral  region 
of  the  NW  and  the  fully  depleted  surface  is  defined  as  (f>m. 
This  convention  should  be  distinguished  from  the  surface 
potential  <Fc  defined  by  <i>s=Ec-Ef  in  a  manner  consistent 
with  the  usage  of  other  authors.  ’  An  approximate  expres¬ 
sion  for  4>m  is  derived  by  assuming  an  abrupt  depletion 
boundary  and  using  Gauss’  law  and  the  cylindrical  symmetry 
to  first  calculate  the  radial  electric  field  in  the  space-charge 
region  and  then  computing  q<fim,  the  work  required  to  move  a 
test  charge  from  rQ  to  Rm.  This  simple  temperature  indepen¬ 
dent  analysis  gives 


.  _qNd 

<f>m~  2s 


Rl-rl 


z  i  | 

-  r„  In  |  — 

rn 


(1) 


In  the  preceding  equation,  q=  1.602  X  10  19  C  and  s  =  sQsn 
where  sQ= 8.854  X  10-14  F/cm  and  sr=  8.9. 
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FIG.  7.  Transcendental  solution  (solid  line)  of  Eq.  (1)  showing  an  example 
of  surface  band  bending  4>m  as  a  function  of  conducting  channel  radius  rQ  (as 
defined  in  Fig.  6)  for  a  presumed  cylindrical  NW  with  Nd=l  X  10“ 16  cm-3 
and  Rm=  150  nm.  The  dashed  line  shows  the  surface  band  bending  cf)p  de¬ 
rived  using  the  planar  approximation  that  may  be  found  in  many  books  (Ref. 
32).  The  figure  illustrates  that  if  4>m  is  known  from  a  separate  measurement 
then  rQ  can  be  calculated  for  a  NW  that  is  not  fully  depleted  provided  the 
constraints  of  the  temperature  independent  model  surrounding  the  discus¬ 
sion  of  Eq.  (1)  are  obeyed. 

We  mention  that  Eq.  (1)  is  a  special  case  of  the  more 
general  (but  still  approximate)  cylindrical  analysis  given  by 
Dobrokhotov  et  al. 14  in  their  work  on  GaN  NW  gas  sensors 
functionalized  with  metal  nanoparticles.  An  example  solution 
of  Eq.  (1)  as  a  function  of  rQ  is  given  in  Fig.  7,  where  rep¬ 
resentative  values  A^=7X1016  cm-3  and  Rm=150  nm  are 
used.  Figure  7  also  illustrates  that  the  band  bending  becomes 
approximately  parabolic  as  the  depletion  layer  becomes  thin 
compared  with  the  NW  diameter  as  expected  from  the  stan¬ 
dard  planar  depletion  model.32 

The  upward  bending  of  the  conduction  band  for  GaN 
surfaces  has  been  estimated  by  electron  photoemission,  scan- 

31  33  34 

ning  Kelvin  probe  microscopy,  and  theoretical  work.  ’  ’ 
From  the  results  of  those  studies  we  expect  that  (fim  should 
fall  somewhere  in  the  range  of  roughly  0. 1-1.1  eV.  Addition¬ 
ally,  Ref.  33  reports  that  there  should  be  no  significant  dif¬ 
ference  in  surface  band  bending  between  polar  c-plane  sur¬ 
faces  and  nonpolar  a- plane  surfaces  for  MOCVD  grown 
GaN  films.  Notably,  the  surface  band  bending  for  GaN  can 
be  strongly  influenced  by  surface  treatments33  and  Ga/N  flux 
ratio  in  the  case  of  MBE  growth.31  This  latter  point  is  sig¬ 
nificant  for  the  present  discussion  since  the  PAMBE  process 
that  enables  the  NW  growth  morphology  requires  a  much 
higher  N  to  Ga  flux  ratio  than  what  is  typically  used  for 
conventional  MBE  growth  of  GaN  films.5  Given  this  range 
of  estimates  reported  in  the  literature,  there  is  no  a  priori 
reason  to  presume  that  (fim  for  our  samples  that  are  examined 
in  ambient  air  should  be  as  high  as  1 . 1  V,  as  reported  in  Ref. 
33,  or  necessarily  ^0.55  V,  as  reported  for  the  c-axis 
PAMBE  grown  GaN  NWs  in  Ref.  13  since  the  latter  are  not 
grown  under  exactly  the  same  conditions  as  ours. 

Given  its  convenience  and  simplicity,  one  is  particularly 
interested  in  the  validity  of  Eq.  (1)  as  <pm^kT/ q,  and  how 
well  the  cylindrical  representation  of  the  hexagonal  NW 
cross  section  applies  in  this  limit.  We  examined  these  issues 
as  follows:  A  temperature  dependent  numerical  Poisson 
solver  applied  to  a  NW  geometry  of  hexagonal  cross  section 


was  used  to  compute  nL  (the  integrated  free  charge  per  unit 
length  of  a  NW)  for  the  representative  cases  including  4>m 
=  0.2,  0.4,  0.6,  0.8  V;  A^=4X1016,  8  X 1016,  and  1.2 
X  1017  cm-3;  and  Dj= 220,  280,  320  nm.  In  these  simula¬ 
tions,  the  acceptor  density  Ns  was  taken  as  uniform  over  the 
NW  surface  and  its  value  was  adjusted  to  compute  the  values 
of  cf)m  just  listed  (at  the  tangential  points  where  Dt  intersects 
the  hexagonal  faces).  The  computed  span  of  Ns  ranged  from 
2X1011  cm-2  (smallest  Db  Nd,  and  (f>m)  up  to  8 
X  1011  cm-2  (largest  Db  Nd,  and  <f)m ).35  Moreover,  given  the 
estimates  of  minority  carrier  lifetime  and  absorption  coeffi¬ 
cient  previously  discussed  in  the  context  of  Figs.  3(c)  and 
3(d)  we  conclude  that  P(325)  in  the  range  of 
102-103  W/cm2  should  be  sufficient  to  nearly  or  fully 
screen  the  depletion  field  over  these  representative  spans  of 
Db  Nd,  and  (f>m.  The  simulation  also  reveals  that  <fim  is  a 
strongly  nonlinear  function  of  Ns  and  for  Ns~2 
X  1010  cm-2,  the  integrated  free  charge  per  unit  length  nL  is 
within  90%-95%  of  the  flatband  value.  Thus  we  expect  that 
a  density  of  photoexcited  holes  nph  of  roughly  2 
X  1010  cm-2  will  yield  — 

The  span  of  values  of  <fim  and  Nd  listed  in  the  previous 
paragraph  was  then  used  as  input  to  calculate  nL  by  means  of 
Eq.  (1).  In  these  calculations  Rm  was  varied  until  the  com¬ 
puted  values  of  nL  equated  with  those  obtained  using  the 
Poisson  solver  for  the  examples  of  Dt  listed.  Dm  calculated  in 
this  fashion  represents  the  diameter  of  an  “equivalent”  cylin¬ 
drical  NW  and  we  label  this  diameter  as  Dem.  We  thus  arrived 
at  the  values  for  Dem  of  Dem=Dt+ 13  ±  5  nm,  Di+10±3  nm, 
and  D/+9±4  nm  that  correspond,  respectively,  to  the  three 
cases  of  Nd  just  listed.  Hence,  over  this  representative  span 
of  <f>m,  Nd,  and  Db  we  find  that  10  nm.  It  is  also 

instructive  to  compare  the  thickness  of  the  depletion  layer 
derived  in  the  foregoing  analysis  to  the  effective  Debye 
length  L,  where  L=^skT/ q2Nd  (assuming  that  only  donors 
are  present).36  Then  another  statement  of  the  applicability 
range  of  Eq.  (1)  is  that  the  depletion  layer  thickness  Deml 2 
-  rQ  should  be  greater  than  roughly  four  to  five  effective  De¬ 
bye  lengths.  Hence,  the  depletion  in  a  NW  of  hexagonal 
cross  section  for  the  diameters  considered  can  be  reasonably 
represented  by  a  cylinder  of  constant  diameter  Dem  (with 
Dem^sDi+ 10  nm),  and  Eq.  (1)  is  an  approximate  representa¬ 
tion  of  the  depletion  in  this  structure  provided  that  0.2  V 
<<£m<0.8  V  and  4  X  1016  cm'3 < Nd<  1.2 X  1017  cm'3. 

2.  Graphical  analysis  of  steady-state  dark  current  and 
illustration  of  uncertainties 

The  numerical  problem  for  analyzing  the  steady- state 
dark  conductivity  in  the  NWs  involves  first  examining 
samples  with  measurable  dark  current  Id  and  then  calculating 
estimates  of  Nd,  /x,  and  <pm  from 

Id  =  NdM7rr20Va/Lg,  (2) 

where  Lg  is  the  electrode  gap  and  Va  is  the  applied  voltage. 
Under  the  assumption  that  NWs  from  a  common  growth 
batch  have  common  transport  properties  and  surface  band 
bending,  one  might  expect  that  measurements  of  Id  from  at 
least  three  separate  NWs  could  allow  for  the  unique  determi¬ 
nation  of  Nd,  /x,  and  <f)m.  In  practice,  however,  this  is  not 


Downloaded  12  Feb  2010  to  132.163.53.130.  Redistribution  subject  to  AIP  license  or  copyright;  see  http://jap.aip.org/jap/copyright.jsp 


034318-7 


Sanford  et  at. 


J.  Appl.  Phys.  107,  034318  (2010) 


FIG.  8.  (a)  Graphs  of  Nd  vs  mobility  computed  from  Eqs.  (1)  and  (2)  for 
three  hypothetical  NWs  with  diameters  Dm  as  listed  in  the  inset,  4>m 
=0.6  V,  and  each  conducting  known  values  of  dark  current  Id  as  given  in 
the  text.  All  three  curves  intersect  at  a  common  point  corresponding  to  the 
input  values  of  Nd  and  /x.  Also  shown  is  the  collapse  of  the  three  curves  to 
a  single  graph  of  Nd  vs  mobility  when  (f)m= 0  V.  (b)  Graphs  of  Nd  vs  mo¬ 
bility  computed  as  in  (a)  but  with  <fim= 0.7  Y  and  the  inclusion  of  an  uncer¬ 
tainty  in  Dm  of  ±10  nm  for  each  NW.  The  shaded  area  illustrates  the  am¬ 
biguity  in  the  common  solution  for  Nd  and  /x  for  the  three  NWs.  Hence,  if 
the  input  values  of  Id  corresponded  to  measured  data,  the  possible  solutions 
for  Nd  and  /x  given  the  assumed  value  of  4>m  and  uncertainty  in  Dm  cover  the 
entire  shaded  area,  (c)  Graphs  of  Nd  vs  mobility  computed  as  in  (a)  but  with 
0m=O.3  V  and  the  inclusion  of  an  uncertainty  in  Dm  of  ±10  nm  for  each 
NW.  The  result  illustrates  that  the  range  of  possible  solutions  of  Nd  and  /x 
depends  strongly  on  (f>m. 

feasible  since  the  uncertainty  in  the  measurements  of  the  NW 
diameters  is  a  significant  limitation  in  this  analysis  as  we 
now  illustrate. 

Consider  three  hypothetical  two-terminal  cylindrical 
NW  structures  with  respective  diameters  D7^2,3  =  250,  300, 
and  380  nm.  We  set  <fim=0.6  V,  A^=1X1017  cm-3,  /z 
=  300  cm2/(V  s),  Lg= 4  /tm,  and  Va  =  1  V  for  each  NW. 
Equations  (1)  and  (2)  may  be  used  to  numerically  solve  for 
the  respective  dark  currents  conducted  in  these  wires,  yield¬ 
ing  Id23  =  0.3868,  1.4655,  and  4.226  /zA.  Figure  8(a)  illus¬ 
trates  the  locus  of  solutions  of  Eq.  (2)  for  Nd  and  /z  corre¬ 
sponding  to  the  input  values  of  Ild  2,3  and  the  simultaneous 
solution  conforming  to  the  initially  assigned  values  of  Nd  and 
/z.  Also  illustrated  in  the  figure  is  the  computed  flatband 
condition  where  cf)m= 0  and  r0-Rm  for  each  NW,  yielding 
solutions  that  are  necessarily  degenerate  for  the  three  cases. 
In  Fig.  8(b)  we  use  </>m  =  0.7  V  and  illustrate  the  locus  of 
permissible  solutions  for  Nd  and  /z  when  an  uncertainty  in 
diameter  SDm=  ±10  nm  for  each  NW  is  included.  This  re¬ 
sults  in  a  range  of  possible  simultaneous  solutions  of  equal 
probability  with  Nd  spanning  roughly  from  0.94  to  1.12 
X 1017  cm-3  and  /z  spanning  roughly  from  225  to 
425  cm2/(V  s).  The  possible  solutions  for  Nd  and  /z  also 
depend  strongly  on  (f)m ,  as  illustrated  in  Fig.  8(c),  where  we 
show  an  additional  example  resulting  from  (f)m=03  V.  Thus, 
under  the  assumption  that  the  cylindrical  representation  of  a 
NW  is  accurate,  the  “correct”  values  for  carrier  concentration 
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FIG.  9.  Graphical  solution  for  Nd  vs  mobility  over  a  representative  range  of 
4>m  for  samples  1-4  using  the  respective  measured  values  of  Id  and  Dm  as 
given  in  Table  I.  The  uncertainty  in  Dm  is  accounted  for  by  assigning  the 
maximum  and  minimum  values  in  diameter  to  those  of  cylinders  that  cir¬ 
cumscribe  and  inscribe  a  NW  of  hexagonal  cross  section,  as  explained  in  the 
text. 

and  mobility,  i.e.,  A^=1X1017  cm-3  and  /z 

=  300  cm2/(V  s),  cannot  be  recovered  without  significant  er¬ 
ror  even  when  the  SDm  is  only  a  few  percent  of  the  NW 
diameter.  We  now  apply  these  graphical  methods  to  estimate 
ranges  of  Nd  and  /z  for  samples  described  in  Table  I  that 
displays  nonzero  dark  current,  and  then  check  the  consis¬ 
tency  of  the  results  with  the  observed  photocurrent. 

In  Fig.  9  we  illustrate  the  result  for  samples  1-4  where  a 
representative  span  in  (j)m  from  0.2  to  0.6  V  is  considered. 
Table  I  gives  the  corresponding  values  of  Id,  Va ,  Lg ,  and  Dc 
used  as  input  for  each  sample  in  the  calculation.  We  consider 
the  maximum  uncertainty  in  diameter  by  calculating  results 
using  estimated  minimum  and  maximum  values  of  diameter 
for  each  sample.  The  extrema  in  Dm  are  estimated  as  follows: 
FESEM  was  used  to  measure  Dm  in  plan  view.  The  NWs 
typically  rest  with  a  hexagonal  facet  in  contact  with  the  sub¬ 
strate  (or  bottom  electrode  pad).  Thus,  £>m~Z)c,  and  values 
of  Dc  appearing  in  Table  I  correspond  to  such  estimated 
maxima.  The  estimated  minimum  NW  diameters  D{  corre¬ 
spond  to  inscribed  cylinders  as  defined  earlier.  As  shown  in 
Fig.  9,  when  this  estimated  maximum  uncertainty  in  the  NW 
diameter  is  taken  into  account,  the  measured  values  of  Id  for 
growth  batch  B738  can  be  accommodated  over  the  given 
representative  range  of  (f)m  with  /z  spanning  from  roughly 
100  to  750  cm2/(V  s)  and  the  corresponding  range  of  Nd 
spanning  from  1.8  X  1017  down  to  4.0  X  1016  cm-3.  Note  that 
Table  I  also  includes  sample  5  that  exhibits  the  lowest  non¬ 
zero  Id  of  batch  B738.  However,  since  this  sample  is  rela¬ 
tively  close  to  depletion,  it  is  poorly  simulated  by  Eq.  (1)  and 
its  measured  value  of  Id  does  not  map  onto  Fig.  9  using  the 
procedure  just  described.  The  samples  of  batch  B992  cannot 
be  analyzed  by  this  procedure  since  there  is  only  sample  9 
thick  enough  (Z)c  =  680  nm)  to  permit  nonzero  Id.  The  ex¬ 
perimental  results  obtained  from  batch  B992  suggest  that  Nd 
for  this  batch  should  be  roughly  ten  times  less  than  for  the 
B738  sample  set.  This  point  will  be  discussed  further  in  Sec. 
II  E. 
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FIG.  10.  (Color  online)  (a)  Schematic  showing  the  reduction  in  surface  band 
bending  from  the  equilibrium  dark  value  cf>m  to  <FL  under  steady  UV  exci¬ 
tation.  Upward  arrows  indicate  photogeneration  of  electron-hole  pairs;  the 
downward  arrows  indicate  aggregate  electron-hole  recombination  via  radia¬ 
tive  and  nonradiative  pathways.  Photoexcited  holes  drift  to  the  surface  of  the 
NW  under  the  action  of  the  depletion  field  and  become  trapped,  thereby 
partially  screening  the  ionized  acceptor  states  and  reducing  the  surface  band 
bending.  This  nonequilibrium  situation  is  depicted  schematically  by  separate 
effective  Fermi  levels  for  electrons  and  holes,  as  illustrated.  The  radius  of 
the  neutral  conducting  region  in  the  dark  r0  increases  to  rL  under  UV  illu¬ 
mination,  yielding  increased  conduction  of  the  NW.  (b)  Schematic  showing 
the  transient  recovery  of  the  surface  band  bending  following  cessation  of 
UV  illumination.  Conduction  electrons  surmount  the  barrier  O(f)  by  thermi¬ 
onic  emission  and  recombine  with  trapped  holes  at  the  surface.  As  the  pro¬ 
cess  continues,  <F(f)  eventually  returns  to  the  dark  equilibrium  value  4>m. 


D.  Analysis  and  simulation  of  persistent  photocurrent 
decay 

1.  Development  of  model 


An  approximate  model  for  the  time  decay  of  the  persis¬ 
tent  photocurrent  ID(t)  is  now  developed.  Our  approach  is 

y  37 

motivated  by  the  work  of  Reshchikov  et  al.  in  their  study 
of  UV-induced  transient  photo  voltage  in  planar  GaN  films. 
Under  uniform  above-bandgap  illumination,  photogenerated 
holes  are  swept  to  the  NW  surface  by  the  depletion  field  and 
trapped.  This  accumulation  of  holes  on  the  surface  effec¬ 
tively  compensates  the  negatively  charged  surface  acceptors 
and  reduces  the  upward  bending  of  the  conduction  band  to 
<EL,  where  <f>L  <  <fim.  The  NW  becomes  less  depleted,  and  the 
radius  of  the  neutral  region  increases  to  rL  with  r0<rL 
<Rm.  This  nonequilibrium  situation  is  depicted  in  Fig.  10(a). 
When  the  illumination  is  ceased,  an  instantaneous  initial  per¬ 
sistent  current  I°p  flows  and  the  total  initial  current  in  the  NW 
is  I°p+Id.  Recombination  occurs  via  electrons  at  the  bottom  of 
the  conduction  band  in  the  neutral  region  surmounting  the 
surface  barrier  <EL  by  thermionic  emission  and  recombining 
with  the  holes  trapped  at  the  surface,  as  illustrated  in  Fig. 
10(b).  The  surface  band  bending  increases  with  time  as  <F(f), 
the  radius  of  the  neutral  region  decreases  with  time  as  r0(t), 
Ip(t)  decreases  from  its  initial  value  of  I°p,  and  the  recombi¬ 
nation  rate  slows  because  the  increasing  O(f)  results  in  a 
greater  barrier  to  thermionic  emission.  This  approximate 
model  neglects  other  possible  trap  states,  tunneling  of  carri¬ 
ers  through  the  surface  barrier,  and  additional  radiative  and 
nonradiative  recombination  pathways  that  may  exist.  As  the 
dark  equilibrium  state  is  restored,  <f>(f)— ><£m,  r0(t)-*r0, 
I p(t)  — >  0,  and  the  current  returns  to  Id.  The  rate  of  electron- 
hole  recombination  may  be  expressed  as 


dps(t) 

dt 


=  -  ps(t)CnNd  exp 


(3) 


In  the  preceding  equation  ps  is  the  time  dependent  photoge¬ 
nerated  hole  density  at  the  surface,  Cn  is  the  electron-hole 


(a)  Time  (s)  (b) 


Time  (s) 


FIG.  11.  (a)  Hypothetical  photoconductive  decay  computed  from  Eq.  (4) 
with  Dm= 200  nm,  Vj=3X1017  cm-3,  <fim= 0.65  V,  and  three  values  of 
capture  cross  section  an\  10-16,  10-15,  and  10-14  cm2.  The  result  is  compa¬ 
rable  with  a  similar  treatment  of  transient  surface  photo  voltage  given  in  Ref. 
37.  (b)  Hypothetical  photoconductive  decay  computed  from  Eq.  (4)  with 
Dm= 350  nm,  Nd= 5  X  1016  cm-3,  </>m=0.3  V,  and  three  values  of  an:  10-22, 
10-21,  and  10-20  cm2. 


capture  coefficient,  and  the  density  of  thermally  activated 
electrons  from  the  bottom  of  the  conduction  band  in  the  neu¬ 
tral  region  of  the  NW  that  is  available  for  surface  recombi¬ 
nation  is  given  by  Nd  multiplied  by  the  Boltzmann  factor.  Cn 
is  related  to  the  capture  cross  section  crn  by  the  relation  Cn 
=  crn^8kT/  7 Tmn,  where  mn  is  the  effective  mass  of  conduction 
band  electrons  and  mn  =  0.20me ,  yielding  Cn  =  crn2A 
X  107  cm/s.37 

Equation  (3)  is  readily  put  into  a  form  amenable  to  nu¬ 
merical  simulation  with  the  following  arguments.  From  the 
requirement  of  charge  neutrality  between  the  surface  and  the 
depletion  region,  the  steady- state  density  of  ionized  surface 
acceptors  in  the  dark  may  be  expressed  as  Ns=(Nd/2Rm) 
X(R2-r^).  Similarly,  the  steady- state  population  of  ionized 
surface  acceptors  under  constant  illumination  is  N$ 
=  {Ndl2Rm)iR2m-r2I).  Thus,  the  steady-state  population  of 
photoexcited  holes  under  constant  illumination  is  pL0=Ns 
-Nls.  When  the  illumination  ceases,  the  time  dependent 
population  of  photoexcited  holes  is  then  p(t)=Ns-Ns(t).  The 
time  dependent  undepleted  cross  sectional  area  of  the  NW 
arising  from  photoexcitation  is  defined  as  A(t),  where  A(t) 
=  7r(r2(/)-r2).  Thus  p(t)=NdA(t)  /  (2irRm)  and  Eq.  (3)  is  re¬ 
cast  as 


dAjt) 

dt 


=  -  A{t)CnNd  expl 


(4) 


The  time  dependence  of  the  surface  band  bending  <E(/)  is 
obtained  from  Eq.  (1)  with  rQ  replaced  by  yjr2+A{t)l it. 
Clearly  Ip(t)=NdpLqA(t)Va/ Lg  so  that  the  solution  of  Eq.  (4) 
is  proportional  to  the  time  decay  of  the  persistent  photocur¬ 
rent.  Equation  (4)  is  solved  numerically  using  an  adaptive 
step  fourth-order  Runge-Kutta  algorithm.  The  specific 
Runge-Kutta  coefficients  used  are  given  by  Carnahan  et  al. 

To  put  this  into  context,  we  compute  example  solutions 
of  Eq.  (4)  representing  the  time  decay  of  normalized  persis¬ 
tent  photocurrent  for  an  undepleted  NW  where  Nd,  (f>m9  and 
c rn  conform  to  values  of  these  parameters  used  in  the  analysis 
of  transient  photovoltage  for  planar  GaN  films  presented  in 
Ref.  37.  These  results  are  illustrated  in  Fig.  11(a)  where 
Dm= 200  nm,  Nd=3Xl011  cm-3,  <f>m  =  0.65  V,  and  three 
values  for  capture  cross  section  crn  are  used  (10-16,  10-15,  and 
10“ 14  cm2).  The  model  is  not  expected  to  apply  when  [Rm 
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- r0(t)]<4-5L ;  thus,  we  calculate  that  the  initial  condition 
at  which  Eq.  (4)  should  reasonably  represent  the  decay  of 
normalized  photocurrent  7(f)  in  this  example  starts  at  the 
point  where  7(f=0)  —  0.45  and  [r0(t)I  Rm]2  =  [Ip(t)+Id]/ If 
=  I(t).  Both  our  model  and  that  of  Ref.  37  can  be  viewed  as 
similar;  however,  there  are  some  clear  distinctions.  Within 
the  constraints  previously  discussed,  our  model  solves  for  the 
full  time  dependent  evolution  of  the  surface  band  bending  for 
a  GaN  cylinder  and  the  corresponding  time  dependent  cross 
sectional  area  of  the  neutral  region.  By  contrast,  Ref.  37  finds 
the  time  evolution  for  the  surface  band  bending  for  a  planar 
GaN  film  within  the  approximation  that  the  photoinduced 
variation  in  band  bending  is  small  compared  to  the  dark  sur¬ 
face  band  bending.  Hence,  one  would  expect  similar,  but  not 
directly  comparable  results  for  the  same  input  values  of  Nd 
and  crn  applied  to  both  situations.  In  particular,  in  Fig.  11(a) 
we  find  that  the  persistent  photocurrent  should  decay  in 
roughly  50-100  s  when  crn  =  10-15  cm2  which  is  qualitatively 
consistent  with  the  temporal  evolution  of  surface  band  bend¬ 
ing,  as  calculated  in  Ref.  37.  As  a  second  example,  Fig.  11(b) 
illustrates  the  case  for  Dm= 350  nm,  Nd= 5  X  1016  cm-3,  and 
cf)m=03  V,  and  crn  in  the  range  from  10-22  to  10-20  cm2  are 
used.  The  starting  time  of  the  simulation  was  adjusted  to 
account  for  constraint  imposed  by  L  as  previously  described. 
Of  course  our  key  task  is  to  consider  the  inverse  problem, 

1. e.,  fitting  the  simulation  to  the  photocurrent  decay  data  in 
an  effort  to  solve  for  Nd  and  crn.  The  ideal  uniqueness  of 
these  fitted  solutions  is  discussed  in  the  Appendix. 

Fitting  the  temporal  decay  of  the  persistent  photocurrent 
for  wires  with  nonzero  Id  offers  a  means  to  further  restrict 
the  estimated  acceptable  ranges  of  Nd  and  /z  as  follows.  For 
a  given  NW  with  measured  Id,  I°p,  Ip(t ),  and  Dc,  the  ratio  of 
dark  current  to  Hatband  current  is  given  by  Idl If=(r0/ R{)2. 
From  the  earlier  estimates  of  the  range  of  nL  is  necessary  to 
attain  90%-95%  of  Hatband  conditions,  and  using  Figs.  3(c) 
and  3(d),  If  falls  into  the  approximate  range  of  2.5/ p<  If 
<5 Ip.  The  resulting  span  of  IpHf  is  then  0.2-0.4,  and  we 
may  compute  the  corresponding  bounds  on  r0(t)  given  Id  and 
Dc.  Using  these  quantities  and  Dem,  we  vary  Nd  and  crn  in  a 
least-squares  minimization  procedure  to  fit  Eq.  (4)  to  7(f)  for 
the  case  at  hand. 

2.  Analysis  of  samples  from  batch  B738 

Figures  12(a)-12(d)  illustrate  fitting  the  decay  of  7(f)  for 
samples  1  and  2.  Here  the  data  were  normalized  and  the  start 
times  tQ  for  the  simulations  were  constrained  as  previously 
described.  The  start  time  of  the  simulation  shown  in  Fig. 
12(a)  for  sample  1  (the  thickest  samples  from  batch  B738) 
occurred  at  to  =  0  (at  the  start  of  the  actual  photoconductive 
decay)  under  the  assumption  that  7(f=0)  =  0.2.  In  Fig.  12(b), 
however,  the  fit  of  the  same  data  was  performed  for  7(f=0) 
=  0.4  and  tQ  was  necessarily  delayed  until  0(0  had  recovered 
sufficiently  such  that  the  cylindrical  depletion  model  was 
valid.  Figures  12(c)  and  12(d)  show  similar  fitting  of  7(0  data 
for  sample  2.  Analysis  of  samples  3  and  4  gives  similar  re¬ 
sults  but  utilizes  less  of  the  available  data  since  tQ  must  be 
delayed  further  in  order  to  accommodate  the  limitations  of 


FIG.  12.  (a)  Fit  of  the  simulated  photoconductive  decay  to  the  data  of 
sample  1  is  shown  as  the  solid  line  and  the  normalized  data  are  shown  as 
open  symbols.  The  start  time  of  the  simulation  conforms  to  the  start  of  the 
photoconductive  decay  in  the  approximation  that  7(f=0)  =  0.2.  The  optimum 
solution  is  found  for  Nd= 6.0  X  1016  cm-3  and  crw=2.0X  10-19  cm2,  (b)  Fit 
of  the  simulated  photoconductive  decay  to  the  data  of  sample  1  in  the 
approximation  that  7fi=0)  =  0.4.  This  yields  Nd= 4.1X1016  cm-3  and  crn 
=7.0  X  10-23  cm2,  (c)  Fit  of  the  simulated  photoconductive  decay  to  the  data 
of  sample  2  in  the  approximation  that  7(f=0)=0.2,  yielding  Nd=5.4 
X  1016  cm-3  and  cr?J  =  6.0X  10_2°  cm2,  (d)  Fit  of  the  simulated  photocon¬ 
ductive  decay  to  the  data  of  sample  2  in  the  approximation  that  l(t=0) 
=  0.4,  yielding  Nd= 3.2X  1016  cm-3  and  crn= 3. OX  10-23  cm2. 

the  depletion  approximation  imposed  by  these  progressively 
thinner  samples.  Sample  5  is  too  thin  to  apply  the  model 
even  though  its  Id  A  0. 

The  simulations  reveal  that  the  estimated  uncertainty  in 
setting  the  normalization  of  Ip(t)  results  in  the  fitted  solutions 
for  Nd  and  crn  spanning  roughly  (3-6)  X 1016  cm3  and 
10“23-10“19  cm2,  respectively.  These  results  for  Nd  taken 
together  with  the  (dark)  current-voltage  data  presented  in 
Fig.  9  suggest  that  /z  should  fall  in  the  range  of  roughly 
500-700  cm2/(V  s).  This  result  is  consistent  with  our  prior 
assertion  that  /z  for  UID  MBE  grown  GaN  NWs  should  be 
comparable  to  /z  extrapolated  from  Hall  mobility  data  re¬ 
ported  for  crack- free  MBE  GaN  films  with  comparable 
n-type  background  that  was  grown  on  Si(lll).26 

The  fully  depleted  NWs  from  this  batch  (samples  6-8 
with  Dc  spanning  105-113  nm)  typically  showed  faster  pho¬ 
toconductive  decay  [as  illustrated  in  Fig.  13(b),  sample  8,  for 
the  case  of  no  AFD  coating];  however,  these  samples  are  too 
thin  to  permit  fitting  the  simulation  to  the  data  using  the 
foregoing  procedure.  The  Ip  versus  R(360)  data  illustrated  in 
Fig.  4(b)  for  sample  7  is  generally  representative  of  these 
depleted  NWs  from  this  batch  and  we  again  expect  that  Ip  is 
within  90%-95%  of  If  for  R(325)  =  10  W/cm2  and  that 
0.2If<Ip<0AIf  for  R(360)  =  3.6  W/cm2. 

According  to  Ref.  37,  our  estimated  range  of  crn 
(10“19-10“23  cm2)  should  be  comparatively  low  for  semi- 
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FIG.  13.  (Color)  (a)  Comparative  PPC  of  sample  1  before  (“Bare”)  and  after 
(“ALD”)  coating  with  roughly  10  nm  of  ALD  deposited  A1203.  The  corre¬ 
sponding  excitation  wavelengths  are  indicated,  (b)  Comparative  PPC  of 
sample  8  before  (“Bare”)  and  after  (“ALD”)  coating  with  roughly  10  nm  of 
A1203.  The  corresponding  excitation  wavelengths  are  indicated. 


conductors;  however,  even  smaller  values  (~10-25  cm2) 
have  been  reported  in  pulsed  voltage  studies  of  Si  surfaces 
subjected  to  various  etchants.39  Moreover,  Hirsch  et  al .40  in 
their  study  of  photoconductive  decay  in  ~3  jam  thick  GaN 
films  with  Nd=  1  X  1017  cm-3  reported  room  temperature  de¬ 
cay  times  of  roughly  103  s,  inferred  a  capture  barrier  of 
^0.2  eV,  and  remarked  that  the  capture  cross  section  must  be 
“extremely  small.”  If  we  interpret  their  room  temperature 
data  as  arising  from  the  mechanism  whereby  the  variation  in 
surface  depletion  thickness  accounts  for  the  photoinduced 
few  percent  variation  in  conductivity,  then  given  a  decay 
time  of  103  s  and  associating  their  capture  barrier  with  our 
c t>m ,  we  estimate  that  their  capture  cross  section  is  on  the 
order  of  10-26  cm2.  The  variation  in  conductivity  they  re¬ 
ported  is  not  unexpected  using  our  model  since  the  depletion 
layer  thickness  in  the  dark  under  these  conditions  would  be 
roughly  44  nm,  i.e.,  within  a  few  percent  of  the  film  thick¬ 
ness. 

E.  Effects  of  ALD  coating  on  NW  conductivity  and 
examples  of  below-gap  PC  behavior 

Figures  13(a)  and  13(b)  illustrate  the  comparative  pho¬ 
toconductive  behavior  under  illumination  at  360,  380,  400, 
and  550  nm  for  samples  1  and  8  before  and  after  these  de¬ 
vices  were  conformally  coated  with  roughly  10  nm  of  A1203 
by  ALD.  ALD  coating  was  performed  using  trimethylalumi- 
num  and  H20  as  the  reactants.41-43  There  are  significant  dis¬ 
tinctions  to  note  in  comparing  PC  and  PPC  before  and  after 
ALD  coating.  For  example,  prior  to  A1203  coating,  we  gen¬ 
erally  found  relatively  weak  subgap  PC  which  would  be- 
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FIG.  14.  (Color)  Comparative  PPC  of  ALD  A1203  coated  sample  6.  Initial  4 
h  PPC  scan  following  UV  excitation  was  recorded  1  month  following  ALD 
coating.  Follow-up  PPC  scan  under  ambient  light  and  UV  excitation  was 
recorded  17  months  later.  The  sample  was  stored  in  ambient  light  and  not 
exposed  to  UV  between  the  two  scans. 


come  more  prominent  after  ALD  coating.  Additionally,  the 
decay  rate  of  the  persistent  current  for  a  given  device  prior  to 
coating  was  comparable  for  all  illumination  wavelengths  that 
yielded  PC,  even  though,  as  seen  in  Figs.  13(a)  and  13(b), 
there  was  a  considerably  reduced  photoresponse  for  the  ex¬ 
citation  wavelengths  greater  than  360  nm.  However,  after 
A1203  coating,  the  photoconductive  response  at  all  excitation 
wavelengths  and  the  decay  time  were  substantially  increased. 

An  example  of  the  long-term  aging  behavior  of  PPC  for 
sample  6  coated  with  A1203  is  illustrated  in  Fig.  14.  Prior  to 
ALD  coating,  this  sample  showed  no  measurable  photocur¬ 
rent  under  ambient  light  illumination.  The  photoconductive 
decay  time  for  this  sample  prior  to  ALD  coating  was  com¬ 
parable  to  that  displayed  by  sample  8  in  Fig.  13(b).  UV  PC 
experiments  were  performed  on  sample  6  roughly  1  month 
after  the  ALD  coating.  It  was  then  stored  for  17  months  in 
ambient  light  conditions.  Figure  14  illustrates  the  UV  photo¬ 
response  and  PPC  decay  recorded  within  1  month  of  ALD 
coating  and  the  photoconductive  behavior  taken  after  17 
months.  The  latter  starts  with  the  photocurrent  recorded  un¬ 
der  ambient  light,  as  shown,  followed  by  exposure  to  UV 
light.  When  the  UV  light  was  switched  off,  the  PPC  decay  in 
the  dark  was  initially  much  faster  than  that  illustrated  in  the 
initial  experiment.  However,  the  photocurrent  measured  un¬ 
der  ambient  light  conditions  is  not  recovered  after  7  days  of 
monitoring  Ip(t)  in  the  dark.  Figures  15(a)  and  15(b)  illus¬ 
trate  comparative  PPC  under  360  nm  exposure  before  and 
after  ALD  A1203  coating  of  samples  9  and  10  from  growth 
batch  B992;  they  show  that  the  A1203  coatings  produces  a 
comparable  change  in  behavior  of  photoconductive  response 
for  NWs  from  both  growth  runs. 

Table  I  tabulates  data  for  Ip  [photocurrent  for  P(360) 
=  3.6  W/cm2  before  ALD  coating]  with  /ALD  [photocurrent 
for  P( 360)  =  3.6  W/cm2  after  ALD  coating]  for  samples 
from  batch  B738.  Combining  these  results  with  the  previous 
estimates  of  If  for  these  samples  reveals  that  the  effect  of 
ALD  A1203  coating  increases  the  photocurrent  such  that 
7ALd~^/  for  P(360)  =  3.6  W/cm2.  ALD  A1203  has  been  re¬ 
ported  to  introduce  interface  trap  states  for  GaN  whose  den¬ 
sity  depends  on  processing  and  annealing  conditions.44  Our 
results  showing  increased  conductivity  and  a  slowing  of  the 
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FIG.  15.  (a)  Comparative  PPC  under  UV  illumination  before  and  after  ALD 
A1203  coating  of  sample  9.  (b).  Comparative  PPC  under  UV  illumination 
before  and  after  ALD  A1203  coating  of  sample  10. 


decay  rate  of  PPC  in  the  A1203 -coated  NW  devices  suggest 
that  the  oxide  layer  introduces  long-lived  hole  traps  such  that 
the  surface  band  bending  is  substantially  reduced  under  UV 
illumination  and  recovers  far  more  slowly  (if  at  all)  in  the 
dark  compared  to  the  uncoated  cases.  By  contrast,  studies  of 
comparative  PPC  on  GaN  NWs  under  vacuum  and  exposed 
to  air  suggest  that  adsorbed  oxygen  plays  a  role  in  acceler¬ 
ating  PPC  decay.45 

We  generally  observe  consistent  conductivity  and  PC  be¬ 
havior  for  NWs  taken  from  a  single  growth  run  but  differ¬ 
ences  in  behavior  between  growth  runs.  On  the  other  hand, 
the  low  strain  and  defect- free  nature  of  the  NWs  is  main¬ 
tained  from  run  to  run.  Such  trends  were  evident  for  the  two 
UID  ft -type  sample  batches  considered  in  this  paper  and  for 
the  study  of  Si-doped  ft-type  NWs  published  separately.  We 
would  expect  differences  in  UID  ft-type  behavior  between 
sample  batches  because  it  is  difficult  to  duplicate  exactly  the 
same  PAMBE  growth  conditions  between  growth  runs.  Fur¬ 
thermore,  one  would  expect  that  such  differences  would  be¬ 
come  more  evident  for  samples  with  relatively  low  back¬ 
ground  carrier  concentration.  Drawing  from  our  earlier 
arguments  based  on  the  high  structural  quality  of  these  NWs, 
fi  for  batch  B992  should  also  be  in  the  range  of 
500-700  cm2/(V  s).  Furthermore,  when  these  wires  were 
AFD  coated  with  A1203  we  would  also  expect  that  IALD 
~  If  for  P(360)  =  3.6  W/ cm2.  Taken  together,  these  estimates 
are  consistent  with  Nd~5  X  1014—  1  X  1015  for  this  wire 
batch. 

Other  AFD  coating  schemes  were  also  explored  on  both 
sample  batches.  We  performed  comparative  360  nm  PPC  ex¬ 
periments  on  NW  devices  (samples  11  and  13-16)  before 
and  after  they  were  coated  with  roughly  7  nm  of  TaN  by 
AFD.  The  coating  was  performed  using  ( tert - 
butylimido)tris(diethylamido)  tantalum  and  hydrazine  as  the 
reactants.46  The  results  are  summarized  in  Table  I.  Figures 
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FIG.  16.  (a)  Comparative  PPC  under  UV  illumination  before  and  after  ALD 
TaN  coating  of  sample  15.  (b)  Comparative  PPC  under  UV  illumination 
before  and  after  ALD  TaN  coating  of  sample  16. 

16(a)  and  16(b)  illustrate  the  respective  results  of  PPC  decay 
for  samples  15  and  16.  Since  the  intent  of  these  TaN  coating 
experiments  was  to  obtain  a  qualitative  comparison  with  the 
A1203  AFD-coated  samples,  we  did  not  take  extra  steps  to 
ensure  that  all  test  structures  supported  only  a  single  NW 
between  the  electrodes.  Figure  16(a)  shows  the  PPC  behavior 
for  sample  15  consisting  of  two  NWs  from  batch  B738  (with 
approximate  diameters  of  76  and  160  nm)  that  bridged  a 
common  electrode  gap  of  3.4  jmm  with  —  1  V  dc  bias  ap¬ 
plied.  The  observed  photocurrent  prior  to  AFD  coating  was 
consistent  with  only  the  thinner  wire  conducting  under  UV 
illumination.  The  TaN  AFD-coated  sample  shows  increased 
PPC  compared  to  the  bare  sample,  but  the  effect  is  far  less 
than  that  illustrated  in  Fig.  13(b)  for  Al203-AFD-coated 
samples  of  comparable  diameter  from  the  same  growth 
batch.  Figure  16(b)  illustrates  the  results  for  sample  16  con¬ 
sisting  of  two  bridging  wires  (with  approximate  diameters  of 
150  and  200  nm)  taken  from  batch  B992,  where  the  electrode 
gap  was  4.5  jum  and  +1  V  dc  bias  was  applied.  In  this  case 
we  found  that  the  effect  on  the  PPC  of  the  TaN  AFD  coating 
was  essentially  negligible.  As  noted  in  Table  I,  samples  11, 
13,  and  14  all  showed  an  increase  in  peak  PC  after  AFD  TaN 
coating,  but  we  did  not  observe  any  significant  change  in  the 
PPC  decay  time  compared  to  that  for  the  uncoated  state  of 
these  samples. 

Finally,  we  remark  that  control  experiments  were  per¬ 
formed  to  check  for  possible  conductivity  and  PC  in  both  the 
A1203  AFD  and  TaN  AFD  films.  In  these  experiments,  the 
electrical  continuity  before  and  after  AFD  coating  was 
checked  between  contact  pads  that  were  devoid  of  any  con¬ 
necting  NWs.  In  no  case  was  conductivity  or  PC  of  the  AFD 
deposited  films  resolved  within  the  sensitivity  of  our  mea¬ 
surements. 

III.  DISCUSSION 

Measurement  of  transport  properties  in  semiconductor 
NWs  often  involves  the  fabrication  of  FET  devices  and  esti- 
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mation  of  carrier  concentration  and  mobility  from  the  con¬ 
ductance  and  pinch-off  conditions.47  Some  authors  have  cau¬ 
tioned  that  the  results  of  such  analyses  will  depend  strongly 
on  the  sophistication  of  the  simulation  tools  and  the  model 
employed.  Our  own  observations  of  drift  and  hysteresis  in 
NW  FET  devices  further  illustrate  that  the  application  of 
sophisticated  simulation  tools  does  not  necessarily  guarantee 
reliable  extraction  of  transport  properties.16,20,49  The  methods 
presented  herein  offer  an  alternative  approach  for  estimating 
Nd ,  /z,  and  <fim  and  require  the  fabrication  of  relatively  simple 
two-terminal  NW  devices  rather  than  three-terminal  FET 
structures.  The  combination  of  the  I-V  and  time  decay  mea¬ 
surements  offers  a  useful  “cross-check”  since  the  latter  do 
not  depend  on  mobility.  However,  for  the  dark  I-V  measure¬ 
ments,  one  must  employ  a  range  of  NW  diameters  such  that 
the  set  of  samples  spans  from  fully  depleted  in  the  dark  to 
substantially  undepleted  while  at  the  same  time  conforming 
to  the  limitations  in  the  depletion  approximation  as  described 
earlier.  Of  course,  abandoning  the  convenience  and  simplic¬ 
ity  offered  by  the  cylindrical  depletion  model  of  Eq.  (1)  and 
employing  a  fully  temperature  dependent  numerical  solution 
scheme  should  improve  results.  Nonetheless,  selection  of  a 
sufficient  span  of  Dc  is  still  required  and  could  prove  chal¬ 
lenging  depending  on  the  doping  level  and  growth  condi¬ 
tions.  Indeed,  as  described  above,  for  a  doping  level  in  the 
range  of  1014-1015  cm-3,  it  was  difficult  to  find  NWs  thick 
enough  to  yield  nonzero  dark  current.  On  the  other  extreme, 
a  n-type  GaN  NW  doped  at  ~1018  cm-3  with  <pm=03  V 
would  be  completely  depleted  for  Dm~  50  nm,  and  the 
growth  conditions  may  not  yield  a  sufficient  distribution  of 
diameters  in  this  range. 

The  simulations  of  photoconductive  decay  illustrate  that 
relatively  small  uncertainties  in  Nd  (factors  of  2)  can  result  in 
much  larger  uncertainties  in  crn.  Consequently,  it  is  difficult 
to  directly  compare  our  estimated  results  for  crn  with  those  of 
Ref.  37  even  though  the  formal  analysis  is  quite  similar. 
Additionally,  the  uncertainty  encountered  in  setting  the  nor¬ 
malization  of  the  photocurrent  decay  compounds  these  prob¬ 
lems,  as  illustrated  in  Figs.  12(a)-12(d).  Therefore,  the  theo¬ 
retically  obtainable  fits  that  should  return  unique  values  of 
Nd  and  an ,  as  described  in  the  Appendix,  are  not  obtained  in 
practice;  however,  one  can  use  this  analysis  to  place  bounds 
on  these  quantities.  In  this  manner  we  were  able  to  estimate 
that  3  X  1016  cm_3< Nd<6  X  1016  cm-3  for  batch  B738  and 
then  use  the  data  of  Fig.  9  to  determine  that  /z  should  fall 
into  the  range  of  500-700  cm2/(V  s)  as  one  might  antici¬ 
pate  for  low  defect  density  MBE  grown  GaN.26 

We  acknowledge  the  possibility  that  variations  in  carrier 
concentration  between  wires  or  within  wires  from  a  single 
growth  run  could  certainly  exist.  However,  the  general  trends 
illustrating  the  variation  in  conductivity  with  diameter  (both 
herein  and  in  Ref.  20)  suggest  that  this  is  a  secondary  effect 
if  it  occurs. 

The  subgap  PC  was  observed  on  samples  before  and 
after  ALD  coating  with  A1203.  The  intensity  illuminating  the 
samples  at  these  wavelengths  (380,  400,  and  550  nm)  was 
within  a  factor  of  2-3  of  that  available  at  360  nm.  Several 
minutes  were  typically  required  for  the  photocurrent  to  reach 
equilibrium  under  constant  illumination  at  these  wave- 


FIG.  17.  3D  map  of  the  least-squared  solution  of  Eq.  (4)  using  as  input  data 
the  hypothetical  case  illustrated  in  Fig.  11(b)  with  Dm= 350  nm,  Nd= 5 
X  1016  cm-3,  and  an=  1  X  10-21  cm2.  The  map  is  scaled  such  that  Nd  ranges 
from  4.5  X  1016  to  5.5  X  1016  cm-3  and  crn  ranges  from  5  X  10-22  to  1.48 
X  10-21  cm2.  The  vertical  axis  is  the  negative  logarithm  of  the  least-squared 
deviation  between  the  input  and  computed  normalized  photoconductive  de¬ 
cay  and  is  arbitrarily  truncated  at  a  value  of  30.  The  pole  conforms  to  a 
solution  coincident  with  the  values  of  Nd  and  crn  used  to  generate  the  input 
data.  The  locus  of  minima  appearing  as  a  continuous  “ridge”  traversing 
diagonally  across  the  map  conforms  to  the  solution  of  the  limiting  case 
where  <F(f)^0m. 

lengths.  By  contrast,  PC  under  360  nm  illumination  would 
generally  rise  to  a  maximum  level  within  a  few  seconds  for 
uncoated  NWs.  Figures  13(a)  and  13(b)  illustrate  that  the 
decay  time  of  the  PPC  was  comparable  for  all  excitation 
wavelengths  prior  to  AFD  coating.  This  suggests  that  the 
same  decay  process  (e.g.,  photogenerated  holes  trapped  on 
the  surface  recombining  with  conduction  electrons  thermi- 
onically  emitted  from  the  neutral  region)  is  responsible  in  all 
cases.  Armstrong  et  al.  argued  that  subgap  PC  should  be 
dominated  by  majority  carrier  photoemission  from  a  deep 
level  within  the  depletion  region.  Here,  the  photogenerated 
electron  emitted  into  the  conduction  band  would  be  pulled 
into  the  neutral  region  and  the  remaining  positive  charge 
would  partially  screen  the  depletion  field  since  it  is  in  the 
immediate  vicinity  of  the  surface. 

If  this  mechanism  were  dominant  we  would  not  neces¬ 
sarily  expect  the  same  decay  time  as  observed  for  360  nm 
excitation  nor  would  we  expect  to  see  an  increase  in  decay 
time  following  AFD  coating.  However,  uncoated  NWs  show 
roughly  the  same  PPC  decay  times  when  excited  above  or 
below  gap  and  PPC  decay  times  always  increase  after  AFD 
A1203  coating  for  all  excitation  wavelengths  considered. 
These  observations  suggest  that  the  surface  plays  the  domi¬ 
nant  role  and  that  subgap  excitation  produces  mobile  holes 
via  transitions  to  intermediate  trap  and/or  defect  levels.  This 
process  would  first  involve  majority  carrier  photoemission  of 
a  deep  level,  as  suggested  in  Ref.  18,  followed  by  photoex¬ 
citation  of  an  electron  from  the  valence  band  to  the  same 
deep  level,  thus  leaving  a  mobile  hole  available  to  be  swept 
to  the  surface  and  participate  in  screening  the  depletion  field. 

While  such  defect/trap  states  plausibly  explain  the  ob¬ 
served  subgap  PC  within  the  context  of  our  model,  we  saw 
no  evidence  of  midgap  “yellow  band”  PF  from  these  NWs.50 
PF  was  performed  at  16  and  300  K  on  samples  2,  4,  6,  and  8 
(prior  to  A1203  AFD  coating)  and  on  other  various  single 
NW  samples  from  both  growth  runs  before  and  after  AFD 
alumina  coating.8  In  no  cases  did  we  observe  below-gap  PF 
emission.  Finally,  we  also  remark  that  we  do  not  observe 
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Nd  (scaled) 

FIG.  18.  Plan  view  shaded  relief  of  the  3D  map  of  Fig.  17.  As  explained  in 
the  text,  the  superimposed  dashed  line  is  a  separate  solution  of  the  limiting 
case  represented  by  a  solution  of  Eq.  (Al). 

sensitivity  of  the  photocurrent  to  polarization  of  the  excita¬ 
tion  light,  as  has  been  reported  in  other  GaN  NW  photocon¬ 
duction  studies.51 

IV.  CONCLUSIONS 

Combined  metrology  methods  involving  current- voltage 
measurements,  PC,  and  transient  PC  offer  a  means  to  esti¬ 
mate  carrier  concentration,  band  bending,  drift  mobility,  and 
depletion  layer  thickness  in  GaN  NWs.  The  results  illustrate 
that  primary  sources  of  uncertainty  in  these  studies  are  mea¬ 
surements  of  the  diameters  of  the  NW  samples  and  estimates 
of  the  relative  magnitude  of  the  steady- state  photocurrent 
with  respect  to  the  maximum  flatband  current.  The  method¬ 
ology  requires  an  ensemble  of  uniform  individual  NWs  span¬ 
ning  various  diameters  such  that  samples  range  from  ones 
that  are  nearly  depleted  in  the  dark  to  ones  that  are  strongly 
conductive  in  the  dark.  An  additional  constraint  is  that  the 
depletion  layer  thickness  should  be  greater  than  four  to  five 
effective  Debye  lengths  in  order  for  the  simple  temperature 
independent  cylindrical  depletion  model  to  be  valid.  For  NW 
samples  from  a  single  growth  run  that  conformed  to  these 
constraints,  we  found  that  UID  ft-type  free  carrier  concentra¬ 
tion,  surface  band  bending,  and  electron  drift  mobility  fall 
into  the  approximate  respective  ranges  of  (3-6) 
X  1016  cm-3,  0.2-0.3  V,  and  500-700  cm2/(V  s).  A  second 
set  of  more  heavily  depleted  NWs  from  a  separate  growth 
run  with  an  estimated  10-100  times  lower  free  carrier  con¬ 
centration  could  not  be  fully  analyzed  by  these  methods. 
NWs  were  driven  to  apparent  flatband  conditions  by  continu¬ 
ous  325  nm  illumination  at  intensities  near  or  exceeding 
100  W/cm2.  ALD  A1203  coatings  were  found  to  strongly 
increase  the  PC  such  that  nearly  flatband  conditions  could  be 
attained  at  substantially  lower  illumination  intensities.  By 
contrast,  TaN  ALD  coatings  were  found  to  have  less  pertur¬ 
bation  on  the  conductivity  properties. 
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APPENDIX:  DISCUSSION  OF  THE  UNIQUENESS  OF 
THE  PHOTOCONDUCTIVE  DECAY  FITTING 
PROCEDURE 

We  examine  the  problem  of  fitting  photoconductive  de¬ 
cay  data  as  a  method  to  recover  unique  values  of  Nd  and  an. 
This  is  motivated  by  the  concern  that  as  the  decay  ap¬ 
proaches  equilibrium  and  d>(/L)  — » <f>m9  the  process  can  be  rep¬ 
resented  by  a  characteristic  lifetime  rd,  where 

Td=(NdCnyl  zM<t>JkT)  (Al) 

and  </>m  is  proportional  to  Nd  via  Eq.  (1).  In  this  limit  an 
entire  locus  of  Nd  and  Cn  (and  an)  exists  as  permissible  so¬ 
lutions  for  a  single  measured  value  of  rd.  This  suggests  the 
need  to  fit  all  of  the  available  decay  data  such  that  O(r) 
<4>m  subject  to  the  constraint  [Rm-r0(t0)]  is  greater  than 
roughly  4-5L  in  an  effort  to  escape  ambiguous  solutions. 

As  an  example,  we  use  as  input  data  the  hypothetical 
case  illustrated  in  Fig.  11(b)  where  A^=5X1016  cm-3  and 
an=  IX  10-21  cm2  were  used  to  compute  the  decay  of  the 
normalized  photocurrent.  We  then  fit  these  data  using  a  least- 
squares  routine  that  searches  for  the  minimum  squared  de¬ 
viation  between  the  data  and  the  solution  of  Eq.  (4)  over  a 
span  of  Nd  and  crn.  The  result  is  illustrated  in  the  respective 
three  dimensional  (3D)  and  shaded  plan  view  relief  maps  of 
Figs.  17  and  18.  The  x  and  y  axes  are  the  respective  scaled 
ranges  of  Nd  and  crn  and  the  z  axis  is  proportional  to  the 
logarithm  of  the  sum-squared  deviation.  The  topography  of 
the  map  reveals  that  the  expected  solution  (. Nd=5 
X  1016  cm-3  and  an=l  X  10-21  cm2)  appears  as  a  pole,  but 
that  there  is  also  a  continuous  remnant  local  minimum  ap¬ 
pearing  as  a  locus  satisfying  Eq.  (Al)  for  the  special  case 
under  consideration  (with  rd= 90.12  s  for  this  hypothetical 
example).  The  dashed  curve  appearing  in  Fig.  18  is  a  sepa¬ 
rate  computed  solution  of  Eq.  (Al)  superimposed  on  the  re¬ 
lief  map  to  illustrate  that  the  remnant  continuous  minimum 
revealed  from  the  least- squares  fitting  procedure  indeed  con¬ 
forms  to  the  limiting  case  represented  by  the  solution  of  Eq. 
(Al).  This  result  also  demonstrates  that  the  efficiency  of  the 
least-squares  fitting  procedure  can  be  substantially  improved 
by  restricting  the  search  grid  of  Nd  and  crn  to  only  follow  a 
band  of  these  parameters  in  the  vicinity  of  the  remnant  solu¬ 
tion. 
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